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Abstract Using deep observations of the Chandra Deep Field South obtained with MIPS 
at 24 /jm, we present our preliminary estimates on the evolution of the infrared 
(IR) luminosity density of the Universe from z=0 to z~l, We find that a pure 
density evolution of the IR luminosity function is clearly excluded by the data. 
The characteristic luminosity Lir* evolves at least by (1+z) 3 ' 5 with lookback 
time, but our monochromatic approach does not allow us to break the degener- 
acy between a pure evolution in luminosity or an evolution in both density and 
luminosity. Our results imply that IR luminous systems (Lir > 10 11 Lq) be- 
come the dominant population contributing to the comoving IR energy density 
beyond z ~ 0.5-0.6. The uncertainties affecting our measurements are largely 
dominated by the poor constraints on the spectral energy distributions that are 
used to translate the observed 24 /im flux into luminosities. 
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Introduction 

Deep infrared and submillimeter observations performed with ISO and 
SCUBA in the late 90s revealed a very strong evolution of IR luminous sys- 
tems with lookback time (e.g., Smail et al. 1997, Elbaz et al. 1999). Since 
these surveys were only sensitive to the brightest IR sources, the quantification 
of this evolution has yet remained strongly debated so far. In particular, the 
relative importance of such IR-bright objects compared to less luminous star- 
bursts, and their respective contributions to the total comoving energy density 
at high redshift is still unclear. Spitzer, the new infrared facility of NASA, is 
now providing a unique opportunity to address this issue in more detail. 
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1. Infrared luminosities of MIPS sources up to z ~ 1 

We performed deep observations of the Chandra Deep Field South with 
MIPS at 24 \im down to a sensitivity limit of 80 /LiJy (5a detection). Cross- 
correlating our data with source catalogs from various optical surveys (i.e., 
VIMOS VLT Deep Survey: Le Fevre et al. 2004; GOODS: Vanzella et al. 
2004; COMBO-17: Wolf et al. 2004; Chandra source follow-up: Szokoly et 
al. 2004), we derived the redshift of 2635 objects detected at 24 ^m and mostly 
located at z < 1.2. We believe that this redshift identification is nearly complete 
up to z ~ 1 . 

Using various libraries of IR luminosity-dependent spectral energy distribu- 
tions (SEDs, e.g., Dale et al. 2001, Chary & Elbaz 2001), we derived the total 
IR luminosities of the MIPS sources from their observed 24 fim flux. Given 
the poor constraint on the true SEDs characterizing these sources at high red- 
shift, the typical uncertainty affecting these estimates could reach a factor of 2 
to 3. More importantly, influence of dust temperature could also add a small 
systematic bias (Chapman et al. 2003) that we have not quantified so far. We 
will explore this effect more thoroughly in a forthcoming work. 

Below z ~ 0.5, we find that the MIPS sources are rather modest emitters at 
infrared wavelengths (median LrR=10 10 L ). At higher redshifts (0.5 <z<l), 
luminous infrared galaxies (LIRGs, 10 11 L < Ltr <10 12 L ) become the dom- 
inant population among the sources detected at 24 /j,m, while we also detect a 
significant number of more modest starbursts and spirals. The most extreme 
sources such as the ultra-luminous IR galaxies (ULIRGs, Lir > 10 12 L Q ) still 
remain quite rare up to z ~ 1.2. 

2. Evolution of the infrared luminosity function 

Based on this sample of 24 /um-selected sources, and using the libraries of 
IR SEDs for computing the /c-corrections, we derived monochromatic lumi- 
nosity functions (LFs) at different IR wavelengths and in various redshift bins 
up to z ~ 1 . Figure 1 shows these luminosity functions plotted at 60 /im and 
compared to the local 60 fim IRAS LF from Takeuchi et al. (2003). A strong 
evolution is clearly noticeable. We find that it can not be described by a uni- 
form increase of the density of the local IR galaxy population ("pure density" 
evolution), but requires a shift of the characteristic luminosity L* by a factor 
of (l+2;) 4 0±0 - 5 . Note that the well-known degeneracy between a pure evolu- 
tion in luminosity or an evolution combining an increase in both density and 
luminosity can not be broken at this stage. 




Figure 1. Luminosity functions estimated at 60 /im with the 1/Vmax formalism in different 
redshift bins between z = and z=l.2 (* symbols). The 3a uncertainties are indicated with 
vertical solid lines. Data points can be fitted by a pure luminosity evolution (L* oc (1 + z) 4 " 2 , 
solid curve) or a combination of luminosity and density evolution (L* oc (1 + z) 3 ' 5 , $* oc 
(1 + z) 1 ' , dashed curve) of the local 60 pim luminosity function (dotted curve). Vertical dashed- 
dotted lines correspond to the 80% completeness limit in each redshift bin. 



3. Star formation history up to z ~ 1 

Star-forming galaxies are believed to be the major component of the 24 /im 
source population (Silva et al. 2004). Assuming the calibration between the 
star formation rate of galaxies and their infrared luminosities (Kennicutt 1998), 
we converted the evolution of the luminosity functions into an history of the 
star formation up to z ~ 1 . Results are shown on Figure 2 where the relative 
contributions of normal starbursts, LIRGs and ULIRGs are also illustrated. We 
find that IR luminous systems (Ltr, > 10 11 L Q ) become the dominant popula- 
tion contributing to the comoving IR energy density beyond z ~ 0.5-0.6 and 
represent ~70% of the star-forming activity at z~ 1. These preliminary re- 
sults will be further developed by Le Floc'h et al. (2005, in prep.). 

References 

Blain, A., Smail, I., Ivison, R., Kneib, J.-R, & Frayer, D. 2002, Physics Reports, 369, 1 1 1 



4 



1 — • — • — i — • — • — • — i — • — • — • — i — • — • — • — i — • — • — • — i — • — • — • — i — r 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Redshift 



Figure 2. The star formation history up to z ~ 1.2 as seen by MIPS 24 /um (thick solid line), 
decomposed into the contributions of normal galaxies (Lir < 10 11 Lq), LIRGs and ULIRGs 
(shaded regions). Uncertainties on these estimates are still significant and they are not reported 
for clarity. Filled squares correspond to the contribution of detected MIPS sources, while 
empty diamonds and the empty square are from Flores et al. 1999 and Thompson et al. 2001, 
respectively. Model predictions from Chary & Elbaz (2001, cross-hatched region), Blain et al. 
(2002, dotted line) and Lagache et al. (2004, dashed line) are also plotted. 
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